function My = Ry(th)
My = [cos(th) 0 sin(th);0 1 0;-sin(th) 0 cos(th)];
end